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TECHNICAL NOTE k2k5 


FLUTTER ANALYSIS OF RECTANGULAR WINGS OF VERY 
LOW ASPECT RATIO 

By Robert W. Fralich and John M. Hedgepeth 


SUMMARY 


A flutter analysis , employing s lender-hody aerodynamic theory and 
thin-plate theory, is made for rectangular wings of very low aspect ratio 
with a constant thickness. The spanwise variation of wing deflection is 
assumed to be given by a parabola, and the chordwise variation is allowed 
complete freedom. The results show the variation of flutter speed and 
mode shape with aspect ratio. Comparisons are made with additional 
results obtained by approximating the chordwise deflection shape by use 
of parabolic or cubic curves. The analysis shows that the cubic approxi- 
mation gives good results for a ratio of chord to semiBpan less than 3 * 


INTRODUCTION 


Flutter analyses of wing and tail surfaces of very low-aspect ratio 
are complicated by the presence of large amounts of chordwise curvature 
in the flutter mode. Indeed, the question of what chord-bending degrees 
of freedom are necessary in order to obtain good results is largely unan- 
swered at present. 

Several studies are available that deal with the effects of chord- 
wise variation of deformations on flutter. (For example, see refs. 1 
to 7* ) In these references various configurations of very low aspect 
ratio were analyzed . In reference 1, complete generality of chordwise 
deflection shape was allowed and the spanwise variation of deflection 
was assumed to be linear. This reference does not shed light on the 
chordwise shape of flutter modes, however, because no flutter condition 
was found. In references 2 to 7 limited account is taken of the chordwise 
deflections either by the superposition of several modes, by using polyno- 
mial expressions, or by the utilization of a number of discrete stations. 

The present paper is concerned with an analysis of the flutter behav- 
ior of the simple rectangular cantilever plate of very low aspect ratio. 
(See fig. 1.) By making use of slender-body aerodynamic theory and by 
assuming a parabolic spanwise deflection shape, the analysis of this 
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configuration can be made -without restricting the chordwise shape of the 
flutter mode. Since flutter is obtained in this case, the complexity 
of the chordwise deflection shape at flutter can be investigated. Results 
are obtained in the form of flutter boundaries, and the chordwise mode 
shapes for two representative cases are illustrated. In addition, the 
errors caused in the flutter boundaries by approximating the chordwise 
deflection shape by use of parabolic or cubic curves are investigated. 


SYMBOLS 


0 ^, 02 , 0^,0^ constants of integration 


c 


chord 


D plate stiffness in bending, Et^yi2(l - ^ 2 ) 

dQ,d^,d 2 ,d^ constants used in equation ( 29 ) 


E 

F(x,t),f(5,t) 

f(s) 


Young's modulus of elasticity 

chordwise deflection shape 

amplitude of chordwise deflection shape 


K 


f lutter-frequency parameter. 



Pl(x,x) 

P(x,y,r) 

p(x,y,T) 


roots of characteristic equation (l8) 
generalized distributed loading (see eq. (l3a)) 

distributed loading per unit area, positive in z-direction 
lift per unit area, positive in z-direction 


q. 


dynamic pressure. 



s 


semispan 


t 


thickness 
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u 

free -stream velocity 

Vi(T) 

generalized aerodynamic leading-edge shear (see eq.. (l^b)) 

V(y,*r) 

leading-edge shear, positive in z-direction 

w(x,y,T) 

deflection, positive in z-direction 

x,y,z 

coordinate system (see fig. l) 

a, 3 ,7 

parameters defined by equations (21) 

e 

mass-ratio parameter, % ^r- 

2 P™^ 

A 

dynamic-pressure parameter, ^ 

M- 

Poisson's ratio 

1 

nondimensional coordinate, x/s 

n 

total potential, energy of system 

p 

free-stream density of fluid 

^El 

density of material 

T 

time 

^(x,y,z,r) 

perturbation velocity potential 

a> 

flutter frequency 

Subscripts : 


I 

imaginary 

R 

real 
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ANALYSIS 


The flutter analysis contained herein is performed in a manner 
similar to that used in the analysis of static aeroelastic divergence 
in reference 8. The configuration treated, a constant-thickness rectan- 
gular plate of very low aspect ratio, is cantilevered from a rigid wall 
and is located in a fluid flow with a free- stream velocity U. (See 
fig. 1.) 

On the basis of the low-aspect-ratio nature of the configuration, 
the following two assumptions are made: First, since the deflection 

shape would vary in a more complicated manner in the chordwise direction 
than in the spanwise direction, a simple spanwise variation is assumed 
and the chordwise variation is allowed to-be arbitrary. Thus, 


w(x,y,r) = y^Cx,-!") 


( 1 ) 


The second assumption is that slender-body aerodynamic theory may be 
used to find the resulting aerodynamic loads when the plate deforms in 
the shape given by equation (l). In this theory, terms in the linearized 
velocity- potential equation that contain streamwise and time derivatives 
are neglected in comparison with terns containing crossflow derivatives. 
This theory is an extension of Jones 1 steady- aerodynamic theory (ref. 9) 

to unsteady aerodynamics as suggested in references 10, 11, and 2. With 
these assumptions the flutter problem is simplified to the extent that 
an exact solution is possible . 


Aerodynamic Forces 

In slender-body aerodynamic theory, the velocity-potential equation 
for linearized flow reduces to Laplace's equation in the crossflow plane: 


b 2 $ + d 2 0 

dy2 



a It should be noted that references 2, 9 , 10, and 11 consider only 
"pointed" wings._ For a rectangular wing it is necessary to treat explic- 
itly the logical consequences of a nonpointed leading edge. The expres- 
sions for loading derived herein, therefore, differ in detail from, for 
example, those given in reference 10. 
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where 0 is the perturbation velocity potential. The assumptions needed 
in deriving this equation are discussed in reference 11. The boundary 
conditions on the upper side of the xy-plane are given as follows: 


For x < 0, 


for 0 = x 5 c, 


0(x,y,+O,T) = 0 


and for 


0(x,y,+ 0,t) = 0 
(x,y,+0,r) = 




(M ^ s ) 
(kl = B ) 


(3a) 


(3b) 


0(x,y,+O,T) = 0(c,y,+O,T) 


(3c) 


In addition, derivatives of the potential at Infinity must be zero. 

The method used in reference 8 is employed herein to calculate the 
velocity potential 0(x,y,+O,r) on the top surface of the plate. Thus, 


0(x,y,+O,T) 



&r\/s 2 

3t 



00 


The lift per unit area in terms of velocity potential is 


p(x,y,r) = 2pU 


80(x,y,+O,-r) + 1 30(x,y,+O,T) 
3x U 3t 


which, upon substitution from equation (4), becomes 


p(x,y,T) 




2 d 2 ? 

U 3x 3 t 


+ 



( 5 ) 
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As a result of applying s lender-body aerodynamic theory to a rec- 
tangular plan form, the velocity potential jumps from a value of zero 
ahead of the leading edge to some finite value at the leading edge. 
Therefore, a concentrated aerodynamic load acts along the leading edge 
and has a magnitude given by 


V(y,x) = 2pU0(O,y,+O,T ) 

which, upon substitution from equation (4), is given by the following 
equation: 


V(y,T ) =--| q 

5 


dF ( 0 , r ) l dF(0,r) 

dx U 8 t 




- y 2 


(6) 


Structural Equilibrium 

The principle of minimum potential energy is used to derive the 
differential equation of structural equilibrium for the function F(x,t) 
in a manner analogous to that used in reference 8 for the corresponding 
divergence problem. Accordingly, the total potential energy, consisting 
of the strain energy of the plate plus the potential energy of the 
external forces, is given by the equation 



where the subscripts x and y denote differentiation with respect to 
x and y, P is the total distributed loading given by 


P 


= P - Pm* 


d^w 

8t2 


(8) 


and V is the shear at the leading edge. (See eq. ( 6 ).) Substituting 
equation (l) into equation ( 7 ) and integrating yield 
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b£ (&f 

5 W/ 


+ l+sF 2 



+ 




&x - 



x,t)F dx - V 1 ( T ) F (°> T ) 


(9) 


where 


Pl(x, T ) 


V3_(t) 



P(x,y,T)y%y 


V(y,T)y%y 


(lOa) 


(lOh) 


It should he noted that the integrated loadings P-j_ and V-|_ are tempo- 
rarily considered to he unrelated to the deflection F for the purposes 
of performing the variation. 

Use of the calculus of variations to minimize the expression for 
potential energy yields the differential equation 


s£ 

5 


^F 

3x^ 


- 2) s 3 + i«F 
5 dx^ 


Pl(x,T) 

D 


(H) 


and the boundary conditions 

£ & + £ ^ 3 ® =. 0 

5 3x2 5 


1 


(x = 0,c) 


f C* - 5^)35 f- - ^ = 0 

5 ax 5 5 ax D 


s5 a^F 2 
5 a x 5 " 3 


5h)b3 I = 0 


(x = 0) 
(x = c) 


( 12 ) 
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Aeroelaetie Solution 

The partial differential equation and its boundary conditions which 
describe the deflection of the plate due to the generalized loads P^(x,t) 

and V;j_(t) are given by equations (ll) and (12). The generalized loads 

are expressed in terms of the deflections by substituting the expressions 
for aerodynamic loads (eqs. ( 5 ) and (6)) and the inertia loads (included 
in eq. (8)) into equations (10). Thus, 


P-.(x,t) = - 


12 


qs 1 - 


V3OO = 


8% 

8x 2 


jt 

12 


2 8% + -A- ] 

U 8x 8t y2 br 2j 


P t 



8F ( 0 , t ) 
8 x 


, 1 8f(o,t) 
U 8t 


s£ 

5 & t 2 


(13a) 


( 13 *) 


If these values of P-^XjT) and V-j_(t) are substituted into the 

differential equation (ll) and into the boundary equations (12), if the 
resulting equations are made nondimens ional by letting 



and if the time dependence is taken to be 


F(|,t) = f(i)e loyr 


the differential equation becomes 


f IV (S) + 4(A - a^f'U) + 4iK \[^ f '(I) + 


12 


20 - K 2 fl + 

12 


f(6) * 0 (i4) 
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and tbs boundary conditiona are 


f"(6) + agf(l) = 0 


f '"(i) + 4(X - a 5 )f *(&) + 2iK f(6) = 0 

f m (g) - ItejE'Cl) = 0 
where 




5rt qs^ 
48 D 


K = CD 



and 


e 


3t pS 
2 p^t 




( 1 = 0 ), 

M) 


( 15 ) 


(l6a) 

( 16 b) 

(l6c) 


It should be noted that real values of ca (or K) infer simple harmonic 
motion and should yield the flutter boundary. 
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Tiie general solution of the differential equation (li+) can he 
written 


+ Cge 


+ C-e 
3 


where are the roots of the characteristic equation 

m 1 *- - k(\ - aijm 2 - kK ]j^- m + 20 - K 2 (l + =0 (l8) 


Substitution of equation (17) into equations (15) yields 
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The condition for a nontrivial solution is obtained by setting the 
determinant of the coefficients equal to zero. When this determinant 
is expanded and simplified, the following equation is obtained: 


(mi - m 2 )(m 5 - 


* Aj4(Ai2 ♦ 


( m l “ m 3)( m 4 “ m 2) 


A 15 (% 2 ♦ + MAji + B3i)e iMI 


( m l “ " * 3 ) 


A 1 ,(a 2 , + B 23)e 1 (w)f + A 32 (A ta + 


= 0 


( 20 ) 


where 

Apq = Aqp = - ij-apa^ - a2^m p 2 + + mq 2 ) - Ij-ajm-pHiq + m p ^nq 2 

Bpq = Bqp = ^a 2 + 2K (n>p + mq) + 4Ampmq 


Equation (20) is expressed in terms of the roots of the characteristic 
equation (l8) which depend on the values of A, K, and e. Express ing 
this dependence explicitly is impossible; therefore, the following pro- 
cedure is used to find a parametric relationship: Examination of the 

characteristic equation shows that the roots must have the form 


mq = 7 + ig 
m 2 = 7 - ip 
m^ = - 7 + ia 
mij. = - 7 - ia 


( 21 ) 
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where 7 is a positive real quantity and 


a 



. k S 7k 
+ 7 II 12 


(22a) 



| / 

7 2 

- 2^A - a- 

0-fl 

iW 

in which 







K 

1- 

_5 - (a 2 - 

A + a 1; 



I 


7 2 (i + 2i) 

5eA 






12 



(22b) 


(25) 


Here, only real values of K/7 are considered and a and p can be 
either real or pure imaginary quantities. Equations (21), (22), and (25) 
show that the characteristic roots can be expressed explicitly in terms 
of A, 7, and e and that K can be found from equation (25). When 
a and p, given by equations (22), are both real, the flutter equa- 
tion (20), upon substitution from equation (21), becomes 


64a,p 




sin 


ill + ?§ sinh an C0Bh §c + ^2 cosh ccc ginh £c\ 
sa s s. p s s / 


= 0 


(24a) 


where 


D 1 = “ ( 2Ai A 2 + + A 2 B]_) 
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E 1 = ( A 3 + ^ 3 ) 2 A 3 - (**-7 2 + a 2 + P 2 )a^J + 

(AIj. - A) 2a 2 p 2 Alj. - (by 2 + a 2 + P^Ajj 

F 1 = ( A 3 + ^ 3 ) 2a 2 f3 2 A^ - (by 2 + a 2 + P 2 )a^J + 
a 2 p 2 ^A^ - A) 2A^ - (by 2 - + a 2 + P 2 )a^J 


d 2 = "( a 1 b 2 " AaBl) 


E 2 = 7 


( A 3 + ^3) 


ct 2 (A4 - a) A ^8-/ a + 8 ai - 8A + p llj + 2 ( r 2 _ x) £ ^gA 


F 2 = -7 


( a 3 + ab 3) 


a( 87 2 - 8a x - 8A + ^ f§) + 2^ - x) 5 0. 


P 2 (A4 - a) 


a( 8 ^ + 8a, -SA^g)-^- *)£/§* 


in which 
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i 2 = -^U-»f-f^+(r 2 -* + §0f 

B i ■ y t 1 ' >*>* + 6" 2 ' ^)[ a ' + , 1 /§*) 

b 2 - f (i - + (r 2 - >)(» - f 1^) 

A 3 - ^ (i - ix)(l - 2n) - ^ (1 - 11)* + y (1 - H>7 2 + 

A 4 - y (l - n) •+ 7 2 


and the expressions for <x, f 3 , and K/j are given by equations ( 22 ) 

and (25). If a is real and p is imaginary in equations (22), the 
flutter equation ( 20 ) has the form 


L64<x(3^ 


2rc oc ^l c F 1 . , <xc ^l c , 

D cos + e cosh — cos + — — sinh — sin 1 + 

n X S B aPl B £ 


3ln 2zc + ^2 


= 0 


(24b) 


where 


Pi 
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If both, a and 3 are Imaginary, equation (20) becomes 


-64 Oq.Pl 


2 /c 


c^c 3-jC F 2 


D-, cos -1— + En cos cos + 

y x S x 6 S 


Ctrj^C 3^ C \ 

sin sin + 

s s J 


l 2-vc ®2 a l c Pi c Fo ctnC 3i c 

i D 0 sin + _£ S i n C os + — cos sin 
\ 2 s acq s s 3 j. 8 s j 


= 0 


(24c) 


where 


a ! = 





Equation (24a) , (24b) , or (24c) can be solved by the following 
trial-and-error process: First, values of A and £ are assigned; 

then, for varying values of 7 , both the real and imaginary parts of 
the proper form of equations (24) are solved for the corresponding values 
of c/s where they exist. Only for particular values of 7 will these 
two solutions yield identical value s of c/s. By finding these partic- 
ular values of 7 , the appropriate values of c/s associated with the 
assigned values of A and e are determined. Some results obtained 
by this and subsequent procedures are presented and discussed in a later 
section. 


Solution for Large Values of c/s 

The results in the region where c/s is large show that a is 
real and 3 i 8 imaginary. In this region, equation (24b), upon consid- 
eration of the order of magnitude of its terms, reduces to 


Eq cos 




3i c . f Eg 
-=- + xl — cos 
s \ a, 


3}C + Fp 

8 Pi 



= 0 


(25) 


This equation can be satisfied only when 


E 1 F 2- 


^2^1 


= 0 


of- 
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which, when expanded, "becomes 




Only the last 


factor of this equation yields a solution: 



K 

7 



+ 


f 2A + £ 



oc^Aij. = 0 


(26) 


The value of c/s corresponding to the solution of equation (26) can 
he found from equation (25) to he 


c ! 

s ^ 


tan"-*- 



(27) 


Equations (26) and (27) are used to find the variation of A with 
c/s for various values of e in the region where c/s is large. That 
is, value 5 of A and e are assigned and 7 is varied until the two 
values of 7 (there are either two or none of them) are found that sat- 
isfy equation (26). Then, the corresponding values of c/s axe given hy 
equation (27). 


Mode Shapes 


The components of normalized mode shape e"^ 1 " which are in 

f(0) 

phase and out of phase with the maximum leading-edge deflection are given 


hy 


[hill 

and 

f(6) 

[Hoi 

‘r 

f(°) 


, which are the real and imaginary parts, respec- 


tively, of the function 
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im-j^ inigi im^| imv| 

f (| ) + @2 e + C^e + C^e 

f(0) C x + Cg + + Cj^ 


(28) 


Here, "the constants Cg, C^, and are found from the boun dar y equa- 

tions (19) in terms of C]_ to be, respectively, 


C 1 

C2 =4 


m .3 


( m 3 - “4)( a 2 - m i 2 ) A 5h e i ^ 1 ' Hn2 ) s + (m k - mi )(a 2 - 
^)A l 4 e- 1 ( m ^ 3 )l + (mi _ ^ 


c,.2l 
3 A 


{ m k - m 2 ) ( a 2 - m l 2 ) A 2i i . e 1 ^ mi ^ S + ( m l - m 4)( a 2 - 
+ (m2 _ mi) ( a2 _ 


C 1 

Ch = — 
^ A 


Q 

(m z - mjJ kg - m,_ 2 j A 2 je ( 1 B + (1115 - kg - 


♦ ( mi - ^(ag - m^A^M 


where 


A = (mg - m^^a 2 - m^Ag^e ^ 1 5 )s + ^ ^ ^ - m^A^e ( 1 ^ T 


(i% - (a. 2 - mg^A^e ^ 1_H312 ) s 


O |CQ 



Approximate Solutions 


& 


In practical wing construction, an exact solution, such as that obtained in the previous 
sections, would not be feasible and same sort of approximation of the chordvise deflection shape 
would be necessary. One method of approximation is to represent the chordvise variation of 
deflection by the first few terms of a power series. The present configuration (fig. l) is used 
to test the accuracy of such an approximation. 

In this approximate analysis, the assumed chordvise deflection shape 


F(x,t) = [dQ + d^f) + d2 gf + a 3 (f) 5 + • • • 


,in>T 


(29) 


is substituted into equation (9)> the equation for potential energy. Minimization of the poten- 
tial energy with respect to d^, d-j , dg, dj . . . and then substitution of the appropriate 

expressions for P 1 (x,t) and V^t) from equations (lj) yield a set of homogeneous simultaneous 

equations. The solution of these equations is obtained by setting the determinant of the coef- 
ficients equal to zero. Cubic approximation of the chordvise variation of deflection results in 
the following complex equation: , 


t> + 12K 


3 «* 






I u + (la + KX.H + 1** yigr 
£b + [ex + f<2 


5 b + T t ‘= + 


F 12 


^ b + - 3t*)J J '+ ^ b + + ^3 - o’ * ^ + 1 ’T 1 


I b «■ aon J lb + ^(2 - ii)4 + ^ b + Ja + ?P<5 - a)j| + S ^ + i ^ C 7 b + jiy- * + B(3 - 

where 


+ 6(3 - an|f + 12^+ i2!C||— 


o (30} 


b = 


20 - Wl + % 


S. 

s 


1 


StfS-tT N1 VOW 
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By a trial-and-error solution, the flutter-frequency parameter K is 
eliminated and the variation of A with c/s is found for a given value 
of e. The equation for parabolic deformation, obtained from equation (50) 
by deleting the last row and last column, is solved analytically for A in 
terms of c/s and e. 


RESULTS AND DISCUSSION 


Results of the flutter analysis, in which chordwise variation of 
deflection is allowed complete freedom, are given by the flutter boundaries 
in figure 2. For convenience, these boundaries are denoted herein as 
"exact." The variation of dynamic-pressure parameter A with c/s is 
given for two values of mass-ratio parameter e. Here, and in the other 
results presented in this section, Poisson's ratio is taken equal to 1/3- 
Comparison with the static aeroelastic-divergence boundary, obtained from 
reference 8, shows that divergence is less critical than flutter. 

Each flutter boundary, for large values of c/s, consists of a series 
of scallops approaching a limiting value of A. The variation of this 
limiting value of A with e is given in figure 3 . The flutter speed is 
less than the divergence speed and seems to approach it asymptotically for 
high values of e. 

Mode shapes, showing the chordwise variation of deflection at flutter, 
are obtained from the analysis and are shown in figure 4 for two values of 
c/s which correspond to the positions indicated by the tick marks on the 
curve for e = 0.01 in figure 2 . Increasing c/s (decreasing the aspect 
ratio) adds more waves to the mode shape. 

Flutter boundaries, obtained from parabolic and cubic approximations 
of the chordwise deflection shape, are compared with the exact flutter 
boundary in figure 5* Both approximations yield good results for values 
of c/s less than with the cubic curve yielding almost exact results. 
However, for higher values of c/s both approximations yield poor results . 
Apparently, in order to analyze the flutter behavior of wings In this 
region, higher order terms in the chordwise deflection shape must be used. 


CONCLUDING REMARKS 


The present flutter analysis of rectangular plates of very low aspect 
ratio indicates that the flutter -mode shape has an increasing number of 
waves in the chordwise direction as the aspect ratio is decreased. Approxi- 
mation of the chordwise deflection shape by use of parabolic or cubic curves 
yields flutter speeds in fair agreement with those of the "exact" theory 
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for a ratio of chord to semispan less than 5 with the cubic appro ximat ion 
giving almost exact results. For lower aspect ratios, higher order 
approximations must be used. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 17, 4958. 
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